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Modeling methods are developed to predict the discharge characteristics of simulated simplified configura- 
tions for plumbing fixtures connected to horizontal drain branch piping. Computations are carried out to illus- 
trate several methods of determining the effect of various loss coefficients of the drain connection, pipe length, 
pipe diameter, and friction factor. Solutions are obtained for the case of a fixture with a constant head (contin- 
uous refill) and a falling head (emptying of a sink). Numerical solution of the non-linear differential equation for 
the falling head case was obtained by the Runge-Kutta method. Discharge characteristics are presented for a 
range of flows and pipe diameter-to-length ratios representative of plumbing installations. The feasibility of 
developing predictive models for hydraulic characteristics of interconnected plumbing fixture and drainage pip- 
ing systems is shown. The variations of efflux rate with the drain pipe diameter, length, and slope obtained from 
the assumed models are similar in trend to the available experimental data. 
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Nomenclature 

A = cross-sectional area of the container 

a= cross-sectional area of the pipe 

d= pipe diameter 

f= pipe friction factor 

g= acceleration due to gravity 

h = instantaneous elevation of the free surface in the 

container above the center of the outlet of the drain 
h = initial value of h 

h e = initial elevation of the free surface in the container 
above the center of the outlet section of initial pipe or 
inlet section of drain pipe (see fig. 3) 
ki = loss coefficient for the tank drain connection 
& fc = loss coefficient for the initial pipe bends and the 

initial pipe-exit pipe coupling 
i? c = loss coefficient for the sudden change in pipe 
diameter at the location where the initial pipe is 
coupled to the drain pipe 
Q i = a i U i = ideal efflux rate, i.e., the efflux rate when all 

head losses are neglected 
Q e = efflux rate = a e ll e 
u = instantaneous flow velocity 
U e = mean velocity of flow in the exit pipe 
U, : = mean ideal velocity of flow in the exit pipe, that is the 

velocity when all head losses are neglected 
p= pressure 
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P a = atmospheric pressure 
s= tan a = slope or pitch of the drain pipe 
a= drain pipe slope angle 

7 = specific weight of water 

0=(i+A: 1 +yzy^ 

4> l = (k l +K b +f l L,/d!) 
<1> 2 =(f.L./d. + k c ) 

Subscript 

1 refers to initial pipe 
e refers to drain pipe 

Metric Conversion 

feet X 0.3048 = meters 



1. Introduction 

The mathematical modeling for determining the time 
dependent capacity characteristics of the drain-waste-vent 
plumbing system requires that the prescribed wastewater 
loads from a variety of fixtures and devices be determinable 
as the initial and boundary conditions. The development of 
models for calculation of time dependent discharge flows 
from plumbed fixtures as the forcing functions required to 
describe the loads on the drainage system is not well estab- 
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lished. The dynamics of the fixture discharge loading are 
mostly derived from experimental data sources that provide 
the combined loss effects as empirical factors. The theoret- 
ical developments are rendered difficult by the complex 
geometric configuration of the connecting elements of the 
fixture-drainage system (e.g. traps, elbows, tees, overflow 
drain, stopper flow area control, and branch drains) and by 
the turbulent flow conditions and secondary swirling. Fur- 
ther difficulties in a complete theoretical development are 
encountered due to the transition from full pipe flow to par- 
tially filled pipe flow downstream of the trap in the horizon- 
tal branch where changes in the size of elbow fittings, pipe 
diameter and slope often occur. 

The dynamic transient characteristics of connected fix- 
tures are generally assumed to be, a priori, independent of 
the downstream connections, drain or stopper geometry, 
elbow, fittings and the branch piping [1, 2]. 1 Pressure varia- 
tions within the drainage system (propagation of down- 
stream conditions typical of incompressible low flow 
phenomena) are neglected in the usual methods for testing 
fixture performance in open air discharge tests at atmos- 
pheric pressure. The pressure variations, above and below 
atmospheric, occur in multi-story buildings as shown by 
Konen [3] and can influence fixture outflow. 

Requirements for sizing the drain system as embodied in 
the plumbing codes in the U.S. are based upon the concepts 
of "fixture units." This traditional approach for water sup- 
ply pipe sizing was extended to drainage systems, since the 
outflow very nearly equals inflow. This method, however, 
does not introduce the transient dynamics of the discharge 
flow phenomena and is applied as a "quasi-steady-state" 
loading for pipe sizing; it does, however, provide the signifi- 
cant input for peak water loads. 

Test data comparing isolated open air fixture discharge 
characteristics with pipe connected fixtures show that the 
discharge curves are influenced by the connections to the 
drainage system; the predictive models should include the 
relevant constraints due to these connections. The model 
adopted in this study is for a horizontal drain carrying 
waste water from a fixture. The fixture discharge charac- 
teristics are found to be modified when losses due to piping 
and connections are introduced in the flow equations. The 
characteristics of interest are peak efflux rate, efflux rate- 
time history, maximum flow velocity in the horizontal drain, 
velocity-time history and the time required to empty the fix- 
ture. The parameters of interest for fixtures are the volume 
of water in the fixture, geometrical shape of the fixture, fix- 
ture to drain connections, drain exit conditions and physical 
characteristics of the drain pipe (i.e. diameter, length, slope, 
and friction factor). 



1 Figures in brackets indicate literature references at the end of this paper. 



This investigation is intended to be a parametric study 
only, varying each characteristic, one at a time; subsequent 
experiments are required to provide specific loss factors. 

The plumbing fixtures are connected to the horizontal 
drains via connecting elements such as traps, downpipes, 
elbows, tees and other pipe fittings. The fixture discharge 
characteristics are influenced by the pressure losses and 
energy dissipation in these connecting elements, and by the 
physical characteristics of drain pipe. Predictive mathemat- 
ical models and experimental validations are required to 
determine the effects of connecting elements and drain pipe 
characteristics on the fixture discharge parameters. The 
discharge characteristics of various plumbing fixtures with 
drain pipes of different physical characteristics (i.e. 
diameter, length and slope) was reported by Hanslin and 
Perrier [4] in an experimental investigation. Classical ap- 
proaches based upon fundamentals and overly simplified 
models concerning the effects of the variables on the efflux 
rate from fixtures are available in basic hydraulic texts [5 
and 6]. However, these studies for predicting the discharge 
characteristics of plumbing fixtures did not utilize repre- 
sentative models for effects of connector variables on the 
discharge characteristics. 

This report illustrates the development of mathematical 
models for predicting the transient and quasi-steady state 
hydraulic discharge characteristics of a fixture coupled to 
the piping of a plumbing systems. This effort is the first step 
in developing the details of the dissipation effects from 
losses experienced in the flow within plumbing fixtures and 
connecting pipes. The next significant effort is required to 
develop a detailed analytical basis for the partially filled 
horizontal drain pipe. The effects of drain pipe parameters 
on the discharge characteristics of an open container, such 
as a lavatory or sink, under different flow conditions are 
presented. The number of variables for each model studied 
have been reduced in order to illustrate the most significant 
effects. For the steady state efflux rate from a constant head 
(continuous refill) container the effects of drain pipe 
parameters, taken one at a time, are shown. Also the solu- 
tion for emptying a fixture based upon the instantaneous 
efflux rate from a falling head container (emptying fixture) 
is presented. 

2. Unsteady flow drain — Falling head 

Consider an open container with the uniform diameter 
drain pipe of length L fitted with a quick opening valve as 
shown in figure 1 to simulate a lavatory branch drainage 
pipe. Initially the container is full and the valve is closed. 
When the valve is rapidly opened, the available head will 
accelerate the fluid, the velocity of flow and rate of efflux 
will increase from zero to a maximum value. After attaining 
the maximum value, the flow velocity will decrease, because 
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QUCX OPENING VALVE 

FIGURE 1. Lavatory-branch drain schematic. 



the elevation of the water surface in the container will be 
decreasing as the water leaves the container. The final out- 
flow from the tank drain into the pipe would experience a 
discontinuity in the water surface area. 

A one-dimensional unsteady flow analysis is adopted that 
leads to the non-linear equation of motion of the liquid dis- 
charged from a simulated lavatory-branch drainage pipe. 
The one-dimensional method assumes a uniform velocity 
across the flow area. Resistance and dissipation such as 
shear at the wall, turbulent mixing and drain vortex forma- 
tion, and shedding losses are introduced through an empir- 
ical quasi-steady state set of constants. The losses are 
assumed to be proportional to the square of the velocity. 

The equation of motion for the one-dimensional flow is: 



L du 



8 



dt 



= h(t)-$ u 2 /2g 



(1) 



where: h = instantaneous elevation of the free surface 
above the center of the outlet section of the 
drain 

u= instantaneous flow velocity in the drain pipe 

<}> = dissipation factor 



The forcing term on the right hand side of eq (1) causing 
acceleration or deceleration of the fluid is dependent upon 
the decreasing head and the dissipation function. In this 
simplified model the losses accounted for are the friction in 
the pipe, and the loss coefficient for the tank to drain con- 
nection [5, 7]. The function is given by: 



$ = 1 + fc, + fL/d 



(2) 



where: k x = loss coefficient for the tank drain connection 
/= friction factor for the pipe taken as a constant 
factor for the range of Reynolds numbers 
L/d— length to diameter ratio 

The first term of eq (2) simply accounts for the decrease in 
total pressure head as a driving force due to the frictionless 



conversion to the dynamic pressure. The ideal frictionless 
fluid case results from the assumption o£k x =f=0. 

The relationship between the velocity u, at the drain of 
cross-sectional area, a, and the head, h{i), in the tank of 
cross-sectional area, A, is obtained from the continuity 
equation, 



dh 
dt 






(3) 



The minus sign accounts for the falling head in the direc- 
tion of the velocity. 

The simultaneous eqs (1) and (3) are combined by assum- 
ing <}> independent of time, differentiating (1), substituting 
(3) for dh/dt and rearranging the terms to yield 



4> 

(Pu/dt 2 + -j- u 



du 



A L 



(4) 



when 



h <h u a/A = 1. 



The initial conditions are: 



t = o, h = h oi u = o, 



du 
dt 



= gho/L 



with the last condition from eq (1). Equation (4) is recog- 
nized as a non-linear equation with constant coefficients 
when <}> is assumed to take on quasi-steady values. The 
dissipation factors are functions of the local time dependent 
flow conditions; the terms of the function are taken as con- 
stant values over the flow regime, i.e. quasi-steady condi- 
tions, in this report. The solution by a numerical method 
based upon finite difference techniques is readily obtained. 

The solution of eq (4) provides the flow velocity as a func- 
tion of time. The numerical integration of the equation was 
obtained from a computer programmed Runge-Kutta 
method. A numerical integration method provides for 
adjustments in the value of a/A, k, and/at various points in 
the stepwise calculation, if it is desired to account for their 
variations. For a container with constant cross-sectional 
area, a/A changes only when the container empties and 
fluid is still present in the vertical section of the drain, i.e. 
when h is less than h x . 

Example discharge velocity variations with time are 
shown in figure 2. In each of the three samples shown the 
acceleration from rest is very nearly constant over the initial 
time to attain the maximum flow. Since the velocity, u, is 
shown to be linear, it can be inferred from eq (3) that the 
height, h, varies with square of the time, t. The predicted 
discharge characteristics resemble the experimental dis- 
charge characteristics for various plumbing fixtures 
published in the literature [8 and 9]. To obtain quantitative 
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FIGURE 2. Discharge characteristics of a container-drainage system. 



agreement between the predicted and experimental results 
it is necessary to determine the values of various loss coeffi- 
cients from the experiments. As anticipated, the smallest 
time of discharge is obtained for the larger drain area ratio 
when holding the other terms constant; the maximum efflux 
velocity occurs with the lowest dissipation values. 



3. Effects of drainage piping 

The dependence of the discharge characteristics upon the 
interconnections of fixtures and devices with the piping 
system and the fittings at drain connections are investi- 
gated. The schematic for the steady state efflux from a con- 



stant head (continuous refill) open container and drainage 
piping system is shown in figure 3. The assumption of con- 
stant head refill condition simplified the problem to a 
steady state analysis. The governing equations do not re- 
quire the development of computer based numerical solu- 
tions of the differential equations. The effects of the 
variables of the piping system and drain connections can 
then be more readily evaluated in a parametric manner. 
The drainage system consists of an initial pipe of constant 
inside diameter d x and length L u the initial pipe is coupled 
to a drain pipe of constant diameter d e , length L e and slope 
angle a. The mean velocity of the flow (U e ) in the drain pipe 
may be obtained by applying the modified Bernoulli equa- 
tion between "0" and ' V as: 



174 




h + L, tin a 



/V////V//I V7777777 



/ 



FIGURE 3. Constant head container drainage system. 

PJy + Z U 2 J2g = 

(5) 
Pe/y + Z e + U 2 j2g + <f> t U 2 l2g + <f> 2 U 2 J2g 

with the conditions at inlet and outlet given by: 

P = P e = P a ; U o =0; Z -Z e = h e -Lj 

where s = sin a and U l the velocity at the exit of the trap, 
i.e. the inlet of the drain pipe. 
The dissipation functions are (j> 1 and 2 where <f) x is the 
function representing the losses at the tank to drain connec- 
tion and in the initial pipe, and 2 is the function represent- 
ing the losses at the step change in diameters and the 
horizontal drain. The equations are: 

*!=(*! + *»+/; V*i) (6a) 

$ 2 =(f,/L,/d e + k c ) (6b) 

The continuity equation requires that for full bore pipe flow 

U x = PU t , 2 = aja x (7a) 

and the flow rate 

Q. = a t U.. (7b) 

The condition of partially fdled pipe flows is not consid- 
ered here. A more extensive analysis is required to establish 
the transition condition for break away from full flow to par- 
tial filled flow. 



The simplification of terms result in the velocity and flow 
rate equations 



Ljt 



U e = U t [(1 + ~W + P 0! + « 



(8a) 



Q e = <?,- (3* [(1 + -^-)/(l + /3< 0, + 2 )] 1/2 (8b) 

where the symbols are defined in the nomenclature. 

Let U r represent the mean velocity of flow and Q r the ef- 
flux rate when the drainage system of the container of fig- 
ure 3-a consists of just the initial pipe as shown in figure 3b. 
The expressions for U r and Q r , are obtained by letting L e 
and 2 equal to zero and j8 = l as: 



and 



U r = £/,/(! + 0,) 1/2 



Q r = U r a x = Q t l(l + 0J 1 ' 2 . 



(9a) 



(9b) 



Equations (5) through (9) show the dependence of U e and 
Q e on the drain pipe variables (i.e., d„ L„ s, and /,). For a 
falling head container where the conditions are slowly vary- 
ing (i.e. area ratio of pipe to container less than 0.1), and 
may be considered as quasi-steady, these equations also 
give the instantaneous flow velocity and instantaneous ef- 
flux rate. The instantaneous flow velocity and efflux rate 
may be found by replacing h with the instantaneous value 
of the height, h, of the free surface in the container above 
the exit center of the cross-section of the initial pipe. 

Using eqs (8a), (9a), (8b) and (9b), the quantities U t and Q e 
may be nondimensionalized as: 



U./U, = [(l +L e s/h )l(\ +0*0! +0 2 )j 



(10a) 



U./U, = [{1 +LJhJ (1 +0J/U +0 4 0, +0 2 )J (10b) 
and 
Q t /Q i = r 3 2 (UJU 1 ) = r 3 2 {(l +L t s/h )ll + /3 4 1 +0 2 )}"> (n a ) 

Q,/Qr=r3 2 (U./Ur) = 

2 ((1 +L t s/h ) (1 +0^/(1 +18*0, + 2 )}"> 



(lib) 



Equations (6) through (11) are applied to examine the ef- 
fects of variations of any of the drain pipe variables on U 
and Q„ as described in the following sections. 
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3. 1 Effects of variations of the exit pipe diameter, d e 

For examining the effects of d e on U e and Q e the following 
simplification can be made without loss of generality: (1) 5 is 
equal to zero; and (2) drain pipe is very smooth and has a 
short length so that the quantity (f e L e /d e ) is negligible in 
comparison with the quantity k c in dissipation function <j> 2 
that is $2 = k c . 

The value of k c is dependent upon the diameter ratio. For 
P less than 1 (i.e. d e < dj the value of k c has been experi- 
mentally determined by several researchers and may be 
found in the literature, typical values are shown in table 1. 
For /3 > 1 the relationship between k c and fi is given by the 
following expression: 



K c = (p 2 -\y,p>i 



(12) 



Variations of the quantities U e /U„ UJU r , Q e IQi and Q./Q, 
due to variations in d e /d u for 3 different representative 
values of <}>„ are presented in graphical form in figures 4 
and 5. 

Table 1 
Values of loss coefficient, K„ for sudden change in pipe diameter, taken 
from reference 7 



P 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

K. 0.37 0.35 0.32 0.27 0.22 0.17 0.10 0.06 0.02 0.00 



In the case of velocity ratios, in figure 4, it is seen that as 
the diameter ratio increases the velocity ratio decreases 
toward zero. As the drain pipe diameter increases for a 
fixed initial pipe diameter the velocity in the drain pipe will 
decrease or as the initial pipe diameter decreases for a fixed 
drain pipe diameter the velocity in the drain pipe will 
decrease. In figure 5 the efflux rate ratio Q e /Q r goes to 
unity after obtaining a local maximum. The ratio QJQ> 
reaches different asymptotic levels in the limit because the 
higher losses result in lower efflux rate as compared with 
the no loss condition required by Q,-. The critical value of 
diameter ratio ji c , at which Q e /Q r is maximum is determined 
by setting 



d(i 



{Qe/Qr) = 



which yields 



(lc = \T2 



(13a) 



(13b) 



The value of QJQt will also have the maximum value at/J c . 

3.2 Effects of drain pipe length 

To examine the effect of exit pipe length on the efflux the 
pipe slope, 5, is assumed equal to zero and d, equal to d e . 
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FIGURE 4. Water velocity ratio versus drain diameter ratio for efflux from a constant head container. 



176 



1.6 



t M : 



1.2 - 



1.0 - 



2 0.8 - 



0.6 - 



0.4 - 



0.2 - 



0.0 



- 




*1 


Q e /Q r 


Qe/Qj 


- 







1-a 


1-b 





la, 1-b 


\ X~V ' 


2-a 


2-b 


_ 




/ ^\ 3 


3-a 


3-b 





2-a x 








- 




^■_^_2Z2 — ^ — ~m — 





^£&. _^^^^^^ 


- 


3, / 


^ ^2-b 








//// 








- 


l 








x 3-b 






- 


// 








- 


1 


I 1 




1 1 1 1 1 



0.0 



1.0 



2.0 3.0 4.0 5.0 6.0 

DIAMETER RATIO d e /<»1 — - 



8.0 



10 



FIGURE 5. Efflux rate ratio versus drain diameter ratio for efflux from a constant head container. 



By substitution and rearranging terms, the efflux rate 3.3 Effects of drain pipe slope 

equation becomes: 

For this case with diameters d, equal to d e , the equations 
Qe/Qr = {1/[1 + f e L e /d{\ + ii)]} ,/2 (14) may be simplified to illustrate the effects of changes in 

pitch. If the drain pipe is short and very smooth, so that 
The variations of Q e IQ r due to variation in drain pipe quantity f e LJd is negligibly small in comparison to the 
length are presented in figure 6 for these representative quantity(l +<}>,), then the equations are further reduced: 
values of <J>„ and for one value oif e (0.01). 

Qe/Qr = Ue/U r = [(1 + Lj/KW + <|>,)] ,/2 (15) 
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Figure 6. Efflux rate ration versus drain length to diameter 
ratio for efflux from a constant head container. 



and 



Qe/Qr = Ue/Ur = [1 + L e s/h ) 1 ' 2 



(16) 



The variations of Q e IQr due to variations in 5, for three dif- 
ferent values o(Le/h„, are shown in figure 7. The increase in 
efflux with increasing pitch is consistent with the effect of 
the larger gravity force component that accelerates the flow. 
For the longer pipe a greater influence is obtained. 

3.4 Effects of drain pipe friction factor 

For examining only the effects of pipe friction on the flow it 
is assumed that s is equal to zero and d x is equal to d t . The 
simplified equations become: 
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Figure 7. Efflux rate ratio versus pipe slope for efflux from a constant 
head container. 
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Efflux rate ratio versus pipe friction factor for efflux from a con- 
container. 



U./U t = Q./Q r = {l+L.s/h }"* 



(17) 



and 



U./U r = Q./Q T = [{l + W(l + ^+f,Ud t )r* (18) 



The variations of QJQ r due to variation inf., for three dif- 
ferent values of $, and two values of LJd„ are shown in 
figure 8. 

It is noted that for a given value of the friction factor the 
ratio Q e /Q r increases with increasing <{>,. The physical 
interpretation is that increased losses at the tank drain exit 
and the initial pipe section without the added length rapidly 
decrease Q r . In terms of actual efflux rate, Q„ the rate of 
flow is actually falling off with increased losses, but at a 
reduced rate as compared to Q r . 



4. Effects of drain variables on the fall of 

water level and time needed to empty a 

container 

The depth of water and time required to empty the con- 
tainer of figure 3 is considered under quasi-steady flow 
conditions analagous to the region of slowly varying velocity 
shown in figure 2. When the free surface in the container is 
at height, h, above the center of the exit cross-section of the 
drain pipe, the instantaneous mean velocity, U e , is given by 
the relationship: 



The continuity eq(3) may be written as: 
-dh = (a e /A)U e dt. 



(19) 



(20) 



To determine the time required to empty the container, eq 
(20) should be solved rigorously by a succession of integra- 
tions, each over a small range of h and U, with approximate 
mean values for the time interval of the <}>, and § 2 values be- 
ing used. This differs from the approach [6] where along a 
streamline of the flow the head varies with time and a solu- 
tion is obtained based upon the exit pipe velocity. However, 
for examining the effects of the parameters on the time re- 
quired to empty the container, the values of ^ and § 2 are 
assumed constant and the equation is integrated in closed 
form following [6] with the initial condition t=0, h = h , 
which yields 

t = 2(A/a e ) [(1 + fi^ t + <h)/2g] 1 ' 2 Wh. - >fh) (21) 

where t is the time required for the liquid level in the con- 
tainer to fall from h e to h. The time required to empty the 
container, t e , is obtained by setting h=h t (see fig. 3). 

The non-dimensionalized value with respect to ^ is given 
by: 

t./tr = (l/CTU+W.+W/U+M 172 

where t r = {A/ a,) {2/ g y'\\fh e - sfh^l+^f' 1 is the time 
required to empty the container when the drainage system 
consists only of the initial pipe. 

The equations indicate that the time required to empty 
the container will increase when <}> 2 is increased and all of 
the other variables are kept constant To examine the effect 
of the variation of d t on the time required to empty the con- 
tainer, the quantity f&Jd. may be assumed negligible in 
comparison with the quantity (1 + k c ). 

Another form of the equation is: 
(h-h l )l{h.-h l ) = 1 + ri 1 t 2 l{h.-h l )-2-nXylhJ{h t -h l ) (23) 
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where 



5. Conclusions 



(r ] = {aJ2A\2g/(\+W x + $ 2 )Y>\ 



(24) 



Equations (22) and (23) are presented graphically for differ- 
ent representative conditions, in figures 9 and 10, respec- 
tively. Higher efflux rates are anticipated (as has been 
shown) with lower dissipation values, and consequently the 
time for discharge is smaller. The maximum discharge 
value for (i = \f2 results in the minimum time value shown 
in figure 9. 
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FIGURE 9. Time to empty sink versus drain diameter ratio. 



It 



^ ^" n« — 




10 K 12 

TIME |t| IN SECONDS — - 

FIGURE 10. Water level ratio versus time for a falling head container. 



Although the predictive models developed in this paper 
are for simple plumbing systems, the examples demonstrate 
the feasibility of developing mathematical models to predict 
the hydraulic characteristics and performance of complex 
plumbing systems. The theoretical discharge characteristics 
obtained by numerical methods of solution of the governing 
differential equation qualitatively resemble the experimen- 
tal discharge characteristics for various fixtures reported by 
Wyly and Hintz [8] and Pink [9]. To obtain quantitative 
agreement between the predicted and experimental results 
it is necessary to determine the values of various loss coeffi- 
cients from the experiments. The results presented showing 
the predicted variations of efflux rate with the diameter, 
length and slope of the drain pipe are similar to trends to 
the experimental results reported by Hanslin and Perrier 
[4]. 
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